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The effects of ferroic distortion and biaxial strain on the band gap and band edges of SrTiOa (STO) 
are calculated using density functional theory and many-body perturbation theory. Anisotropic 
strains are shown to reduce the gap by breaking degeneracies at the band edges. Ferroic distortions 
are shown to widen the gap by allowing new band edge orbital mixings. Compressive biaxial strains 
raise band edge energies, while tensile strains lower them. To reduce the STO gap, one must lower 
the symmetry from cubic while suppressing ferroic distortions. Our calculations indicate that for 
engineered orientation of the growth direction along [111], the STO gap can be controllably and 
considerably reduced at room temperature. 
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Within the perovskite family (ABO3), a remarkable 
range of functional phenomena — ferroelectricity, ferro- 
magnetism, multiferroicity, and superconductivity, for 
example — can be accessed through chemical substitution 
at the A and B sites. [1] In recent years, biaxial strains 
induced by epitaxy have provided another route to tune 
perovskite properties. [2] In a well-known example, both 
theory and experiment have shown that thin films of 
SrTiOs (STO), a prototypical perovskite insulator that 
takes up a centrosymmetric cubic structure {Pm3m) un- 
der standard conditions, become polar and ferroelectric 
at room temperature for sufficiently large biaxial strains 
induced by coherent epitaxial growth. [H, Q In STO, the 
large biaxial strains achievable via epitaxy provide access 
to a rich landscape of ferroic distortions, both polar fer- 
roelectric and nonpolar antiferrodistortive, [sl, [H, Q that 
are known to alter both ferroelectric polarization and the 
low- frequency dielectric response. [H 

In addition to its well-studied dielectric response, STO 
has also long been explored in solar water splitting ap- 
plications due to its favorable fiat-bandpotentials and 
stability against degradation in water. [7| However, its 
indirect optical gap (3.2 eV [sl) is too large for efficient 
light harvesting in the solar spectrum. Prior work has 
focused on chemical substitution as a means for reducing 
perovskite gaps, [9] but with unintended consequences 
for other properties. In this Letter, we use density func- 
tional theory (DFT) and beyond to examine a new route, 
via anisotropic strain and ferroic distortions, for altering 
the band gap and band edges of STO. We show that 
while anisotropic strain in the absence of ferroic distor- 
tions reduces the gap, ferroic distortions increase it. In 
addition, band edge energies, important for electro- and 
photocatalysis, can also be controllably increased or re- 
duced relative to vacuum. The predicted trends can be 
fully understood with symmetry arguments, suggesting 
well-defined routes for tuning the electronic structures 
of perovskite oxides. For example, using intuition devel- 



oped here, we show that strained thin films of STO grown 
along an alternate orientation, [111], have considerably 
reduced gaps relative to bulk at room temperature. 

STO structural and electronic properties are ob- 
tained using DFT within the local density approxima- 
tion (LDA)^ using the VASP package ^ and PAW po- 
tentials. [11| Electrons taken to be valence are 4s^4p^5s^ 
of Sr, 3523/4523^^2 of Ti, and 25^2/ of O. A plane- 
wave cutoff of 500 eV is used throughout. At room tem- 
perature, STO takes up the cubic perovskite structure 
{PmSm) with a five-atom primitive cell, which we treat 
with a F-centered 8x8x8 /c-point mesh. Near 100 K, 
STO is known to undergo a transition to an antiferrodis- 
tortive phase in which Ti-0 octahedra rotate, doubling 
the unit cell and lowering symmetry to M/mcm. [ill For 
such phases requiring super cells, we use proportionally 
fewer /c-points. These parameters result in excellent con- 
vergence and good agreement with previous work. 

The DFT-LDA near-gap band structure of optimized 
cubic STO (a=3.864 A) is shown in Fig. [H The conduc- 
tion band minimum (CBM) lies at F and consists almost 
entirely of Ti 3d states, while the valence band maximum 
(VBM) lies at R and consists almost entirely of O 2p 
states. As expected within DFT-LDA, the calculated in- 
direct gap of 1.80 eV is far narrower than the experimen- 
tal optical gap of 3.2 eV. |8j To obtain reliable band gaps 
and edges, we use many-body perturbation theory within 
the GW approximation [isj, known to give quantitative 
band gaps and band edge energies for select semiconduct- 
ing systems in the absence of strong excitonic effects, [l^ 
Our computed Go Wo gap of cubic STO (see ^) is 3.82 
eV, 0.6 eV larger than experiment. [28| This overestimate 
is consistent with past calculations for STO [HI and a 
closely-related system, Ti02, where deviations are sug- 
gested to arise from excitonic and polaronic effects. [l6| 
Overestimate notwithstanding, our GoWo band edge en- 
ergies at equilibrium are in reasonable agreement with 
experiment. Experiments place the CBM and VBM 2.6 
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FIG. 1: (Left) DFT-LDA band structure of cubic STO. 
(Right) Schematic pictures of the highest occupied and low- 
est unoccupied states of STO as linear combinations of atomic 
orbit als. Ti atoms are shown in white, O atoms in gray. 



and 5.7 eV below the vacuum level, respectively, [17] com- 
pared with our Go Wo calculated values of 2.48 and 6.30 
eV (see 0). 

Following previous calculations, [H, Q we model biaxi- 
ally strained thin films of STO by periodic calculations in 
which in-plane lattice parameters are fixed while the per- 
pendicular axis and atomic positions are allowed to relax. 
Our structures have in-plane strains of up to ±4%, and 
the resulting lattice parameters (Supplemental Informa- 
tion) and structural phase diagram agree well with previ- 
ous work, [sl, [HI, [gI This range of strain is somewhat opti- 
mistic experimentally, as strains approaching ±4% could 
limit the achievable film thickness, making bulk optical 
measurements difficult, and potentially introducing in- 
terfaces and quantum confinement effects. Nonetheless, 
perovskite substrates that would induce these strains are 
available, [Isl and there are no fundamental limitations 
on stabilization of thin films on these substrates. For 
ultrathin films, our results here represent a quantitative 
basehne upon which future studies of quantum confine- 
ment and interface effects may be understood. 

We start with the most common growth orientation, 
[001], and compute electronic gaps and band edge en- 
ergies for biaxial strains up to ±4%; these are shown in 
Fig. [21 initially for the paraelectric case in which the [001] 
axis is relaxed but further distortion of the internal co- 
ordinates is forbidden {PA/mmm symmetry). Interest- 
ingly, biaxial strains tune both CBM and VBM by several 
tenths of an eV, suggesting a route to optimize each band 
edge relative to the potentials of the water splitting half- 
reactions. From Fig.[2]3, compressive biaxial strains raise 
band edge energies relative to vacuum, and tensile biax- 
ial strains lower these energies. For paraelectric STO, the 
DFT-LDA gap narrows under both compressive and ten- 
sile strain. That the gap is non-monotonic about the cu- 
bic structure is expected on symmetry grounds (Fig.[2t). 
Cubic STO has Oh point group symmetry, and both its 
CBM and VBM are threefold-degenerate. Biaxial reshap- 



ing of the unit cell lowers the point group symmetry to 
^4/i5 splitting these degenerate triples into a spread of 
singles and pairs. As a result, the gap is narrower for 
tetragonal paraelectric STO than for cubic STO. 

Notably, our GoWo corrections to Kohn-Sham DFT- 
LDA eigenvalues (Fig.[2^,b, green) remain quite constant 
with changing biaxial strain. (A similar result was re- 
ported for sihcon under isotropic pressure. [l9|) The same 
is true of GoWo corrections to STO structures with fer- 
roic distortions (Supplemental Information). This sup- 
ports the premise that DFT-LDA, while insufficient for 
producing numerically accurate optical gaps, correctly 
tracks their trends, lol, [2o| For the remainder of this pa- 
per, we restrict our focus to DFT-LDA results, at times 
assuming a constant correction of 1.4 eV — the difference 
between DFT-LDA and room-temperature optical gaps 
of equilibrium STO — to provide concrete estimates for 
gaps under biaxial strain. 

We turn now to the issue of how ferroic distortions af- 
fect the STO band gap. Two classes of commonly studied 
ferroic distortions in perovskites are ferroelectric polar 
modes (FE, relative translation of the cation and anion 
sublattices) and antiferrodistortive modes (AFD, collec- 
tive rotation of the Ti-0 octahedra). While neither dis- 
tortion is present in unstrained STO at room temper- 
ature, both are accessible via strain, temperature, and 
applied fields, [sl, 0, [ill To illustrate how these distor- 
tions might affect the band gap, we freeze in displace- 
ment patterns associated with the eigenvectors of low- 
lying IR-active phonon modes along the [001], [110], and 
[111] axes, and examine the DFT-LDA gap as a function 
of amplitude. We freeze in AFD modes by interpolating 
between cubic STO and optimized structures with AFD 
rotations about the [001], [110], and [111] axes. 

Computed gaps of these distorted STO structures are 
shown in Fig. [3l For modest changes in structural en- 
ergy, both FE (Fig. ^) and AFD (Fig. [3b) modes widen 
the gap by as much as 0.2 eV. This widening can be un- 
derstood with symmetry arguments. As the symmetry is 
lowered under both types of distortion, conduction (pri- 
marily Ti 3(i) and valence (primarily O 2p) states are 
allowed to mix and repel one another. Table H] shows the 
irreducible representations of Ti 3(i and O 2p states at 
F and R within the Oh point group. By symmetry, the 
CBM of cubic STO cannot mix with O 2p states at F or 
R. In the presence of FE distortions, inversion symmetry 
is broken, the distinction between even and odd {g and 
u) parity is eliminated, and the CBM is allowed to mix 
with O 2p states. This mixing was previously illustrated 



by Wheeler et al. [2l| for STO in terms of linear combi- 



nations of atomic orbitals. Similarly, in the presence of 
AFD rotations, doubling of the unit cell renders F and R 
equivalent, allowing the CBM to mix with O 2p states. 
Likewise, the VBM of STO can mix with Ti 3(i states 
only in the presence of FE and/or AFD distortions. As 
Ti 3(i and O 2p states mix, they widen the STO band 
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FIG. 2: a) Computed gaps and b) band edges of paraelectric, biaxially reshaped STO {PA/mmm symmetry) under biaxial 
strain perpendicular to [001]. c) Schematic illustration of the trends in band edges as a function of strain. 
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FIG. 3: Computed gaps (top) of STO structures in which a) 
FE polarization modes and b) AFD rotation modes are frozen 
in. Per- formula- unit structural energies relative to cubic STO 
are also shown (bottom). FE amplitudes are expressed as 
average relative translations of Ti and O atoms, while AFD 
amplitudes are expressed as average translations of O atoms. 
AFD rotations refer to [001] = a°a°c" [110] = a"a"c°, and 
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in Glazer notation. [25 



TABLE I: Irreducible representations of Ti 3(i and O 2p states 
at r and R within Oh point group symmetry. The CBM and 
VBM of cubic STO are shown in red. 



Atomic orbitals 


Irreducible representations 


Ti M at r 


t2g, eg 


Ti M at R 


t2g, eg 


2p at r 


tlu, tlu, t2u 


2j9 at R 


tig, t2g, eg, aig 



gap. While the magnitude of the widening depends on 
the geometry of a given distortion, this trend is general. 
Interestingly, AFD rotations can transform the STO gap 
from indirect to direct, as has been previously noted in 
both theory and experiment. [H, [i^ 



Having established how anisotropic strain and ferroic 
distortions impact the STO electronic structure sepa- 
rately, we present in Fig. |4] the computed gaps of struc- 
tures with ferroic distortions under biaxial strain (assum- 
ing a [001] growth direction). The energetic landscape of 
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FIG. 4: a) Computed gaps of STO structures under biaxial 
strain perpendicular to [001], optimized with individual fer- 
roic distortions, b) Gaps of fully optimized, strained STO 
structures, computed at zero temperature. 



ferroic distortions in ground-state STO has been studied 
extensively within DFT-LDA. [i, i, [26| Consistent with 
this work, full relaxation causes FE and AFD to combine 
along [001] under compressive strain, and along [110] un- 
der tensile strain. All ferroic distortions, optimized sep- 
arately (Fig. |4^), widen the gap relative to paraelectric 
STO. Full atomic relaxation (Fig. HJd) widens the gap 
even further. Using biaxial strain and associated ferroic 
distortions, we predict the optical gap can be tuned by 
nearly 20% under ±4% strains. 

As we have shown, the breaking of symmetry in STO 
has two manifestations with competing effects on the 
band gap — biaxial reshaping of the unit cell that narrows 
the gap through broken degeneracies, and ferroic distor- 
tions that widen the gap through new orbital mixings. 
The data in Fig. HJd suggest that, in biaxially strained 
STO at low temperatures, ferroic distortions have the 
dominant effect on the gap. To engineer a narrower STO 
band gap, we must reshape the unit cell while suppressing 
ferroic distortions. One way to suppress ferroic distor- 
tions is to work at room temperature. Using a Landau- 
Ginzburg-Devonshire model with empirical parameters, 
Pertsev et al. computed a phase diagram of STO under 
[001] biaxial strain. [3j] At 300 K, their model predicts 
that STO is paraelectric between 1.8% compressive and 
1.5% tensile strain. While this provides a range of strain 
within which the STO gap can be narrowed, ferroic dis- 
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FIG. 5: Computed gaps of STO using free energy minimized 
structures at 300 K under biaxial strain perpendicular to 
[111]. All structures are paraelectric. 



tortions widen the gap outside the range. 

If, however, STO is grown epitaxially along an axis 
other than [001], the possibility exists for reshaping the 
unit cell without inducing ferroic distortions. Unlike 
growth along [001], which forces in-plane Ti-0 distances 
to be too short or too long for optimal bonding, a growth 
axis of [111] imposes no constraints on Ti-0 distances. 
We therefore expect that the [111] case does not strongly 
enhance ferroic distortions in STO. Using the same em- 
pirical parameters of Ref. [3] , but applying them to films 
grown along [111] at 300 K, confirms that STO is ex- 
pected to remain paraelectric over the entire range within 
±4% strain. We therefore predict that [111] strain at 
room temperature (achievable via epitaxial growth on 
the (111) surfaces of perovskite substrates) could narrow 
the STO band gap significantly more than [001] strain 
(Fig. [5j). Within DFT-LDA, 4% tensile strain perpen- 
dicular to [111] narrows the Kohn-Sham gap by more 
than 0.3 eV, to 1.46 eV. Using our constant correction, 
we therefore estimate that this approach could lower the 
optical gap of STO below 3 eV. 

To conclude, we have examined the extent to which en- 
gineering of the STO crystal structure can tune its elec- 
tronic gap and band edges. While biaxial strain widens 
the STO gap at zero temperature via ferroic distortions, 
it has the ability to narrow the gap significantly at room 
temperature, especially when applied perpendicular to 
the [111] direction. The observed trends can be un- 
derstood on general symmetry grounds, suggesting that 
similar intuition can be productively applied across the 
broad family of perovskites. 
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